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I)  INTRODUCTION 

The  attempt  to  produce  coherent  millimeter  and  sub^ni  Dimeter  electro¬ 
magnetic  radiation  at  high  power  levels  by  the  deceleration  of  charged  par¬ 
ticles  injected  into  strong  magnetic  fields  is  beset  by  a  number  of  difficulties 
related  to  the  source-detection  problem  in  this  region  of  the  spectrum.  The 
theoretical  study  reported  on  here  was  initiated  to  provide  foundation  and 
guidance  for  the  concurrent  experimental  work. 

Considering  the  complexity  of  the  theoretical  problem,  substantial 
progress  has  been  made  during  the  contract  period.  Approximate  solutions 
have  been  obtained  to  some  of  the  problems  posed,  while  other  problems 
have  been  delineated  more  sharply.  Several  directions  for  further  work 
have  been  uncovered  which  hold  promise  of  being  fruitful .  The  report 
presented  here  summarizes  the  progress  made  during  the  contract  period. 

II)  STATEMENT  OF  THEORETICAL  PROBLEM 

The  general  theoretical  problem  may  be  stated  as  follows:  to  study 
the  dependence  of  gyromagnetic  radiation  production  in  representative 
physical  systems  on  the  system  parameters  in  order  to  define  the  ranges 
of  those  parameters  which  will  give  optimum  performance.  (  Optimum 
performance  being  usually  understood  to  mean  maximum  power  radiated 
within  a  certain  bandwidth  .  )  During  the  course  of  this  contract  one 
particular  type  of  experimental  system  has  been  used,  hence,  the 
theoretical  model  constructed  to  represent  the  physical  system 
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has  been  fixed  within  corresponding  limits.  Two  aspects  of  the  construc¬ 
tion  of  the  theoretical  model  are  dealt  with  in  the  following  two  sections: 
the  choice  of  a  representative  physical  system,  and  the  simplification 
of  the  mathematics  to  obtain  tractable  equations  which  still  yield  mean¬ 
ingful  solutions. 

A)  REPRESENTATIVE  PHYSICAL  SYSTEM 

Basically,  the  theory  must  describe  the  interaction  between  an 
ensemble  of  charged  particles  and  a  strong  mogn.tic  field. 

The  character  of  the  ensemble  may  be  simplified  by  noting  that 
electromagnetic  theory  dictates  that  the  production  of  gyromagnetic 
radiation  is  primarily  a  single  particle  phenomenon,  although  the  char¬ 
acteristics  of  the  emitted  radiation  are  governed  by  the  type  of  magnetic 
field  as  well  as  by  the  mode  of  the  particle  dynamics  in  the  field . 
Furthermore,  if  the  ratio  of  the  mean  free  path  of  the  particles  in  the 
ensemble  to  the  local  gyromagnetic  radius  is  large,  space  charge  effects 
maybe  ignored.  Thus,  to  a  first  approximation,  the  ensemble  may  be 
taken  to  consist  of  non-interacting  particles  and  the  theoretical  model 
then  involves  the  interaction  of  single  particles  with  a  strong  magnetic 
field.  Note  that  in  such  a  model,  coherence  of  the  emitted  radiation 
depends  largely  on  the  initial  relative  phases  of  the  particles  for  a  given 
magnetic  field. 

The  magnetic  field  used  in  the  experimental  arrangement  was  produced 
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by  a  pulsed  magnet  whose  period  was  long  compared  to  the  time  spent 
by  the  particles  in  the  high  field  region,  hence  the  field  could  be  taken 
to  be  magnetostatic .  The  field  possessed  symmetries  characteristic  of 
fields  produced  by  many  pulsed  magnets.  Such  magnets  are  symmetric 
with  respect  to  rotation  about  an  axis  passed  through  the  center  of  the 
air  gap.  Moreover,  they  concentrate  their  fields  within  certain  regions 
along  the  axis  of  symmetry,  frequently  displaying  a  symmetry  about  an 
"equatorial"  plane  passed  perpendicular  to  the  axis  at  the  middle  of 
the  air  gap.  It  is,  in  fact,  the  spatial  inhomogeneity  of  the  field  thus 
produced  that  creates  an  injection  problem  because  the  particles  must  be 
introduced  into  the  field  at  a  point  remote  from  the  high  field  region  in 
order  to  avoid  interference  of  the  magnet  with  the  particle  source.  In 
such  fields,  the  field  magnitude  is  zero  at  +  °>  along  the  axis  and  gradu¬ 
ally  rises  to  a  local  maximum  at  the  magnet  center. 

Thus  the  representative  physical  system  may  be  taken  to  be  a  group 
of  single  particles  moving  in  an  axially  symmotric,  magnetostatic  field 
possessing  a  local  axial  maximum  at  the  equatorial  plane,  and  whose 
magnitude  is  symmetric  with  respect  to  reflections  in  the  equatorial  plane. 

B)  MATHEMATICAL  MODEL 

The  simplifications  necessary  to  establish  a  tractable  mathematical 
model  are  considerably  more  restrictive  in  nature  than  those  introduced 
above.  The  dynamic  assumptions  introduced  are  three:  (  I  )  only  non- 
relativistic  motion  will  bo  considered:  (  2  )  both  the  relativistic 
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and  non -relativistic  radiation  reaction  forces  will  be  neglected,  and 
(  3  )  the  actual  trajectories  will  only  be  dealt  with  in  the  paraxial 
approximation  .  A  discussion  of  these  assumptions  may  be  found  in  the 
paper  entitled:  "The  Injection  of  High  Velocity  Charged  Particles 
into  Strong  Magnetic  Fields,"  scheduled  for  publication  in  the  May, 

1963  issue  of  the  Journal  of  Applied  Physics.  Certain  of  the  qualitative 
results  do  not  require  assumption  (  3  ) . 

The  field  behavior  of  the  magnet  type  outlined  above  may  be 
adequately  represented  for  the  present  purposes  by  a  magnetic  field 
whose  axial  dependence  is  given  by  the  bell-shaped  curve 

B  (  0,  0,  z  )  =  B  /  1  +  (  z/a  )2  (  1  ) 

o 

where  the  z-axis  is  the  axis  of  symmetry  and  a  denotes  the  half-width 
of  the  curve  at  half-maximum.  This  field  has  the  advantage  that  many 
of  its  properties  are  amenable  to  analytic  treatment. 

The  production  of  gyromagnetic  radiation  is,  of  course,  a  relativistic 
effect  in  essence,  hence  due  caution  must  be  exercised  in  the  interpre¬ 
tation  of  results  and  their  application  to  the  study  of  this  phenomenon  . 
Nevertheless,  it  is  felt  that  the  insight  gained  into  the  theoretical 
problem  with  a  non- relativistic  treatment  can  provide  valuable  guidance 
for  the  prosecution  of  further  studies. 

Only  a  few  of  the  basic  properties  of  non-relativistic  motion  in  an 
axially  symmetric  magnetostatic  field  have  been  worked  out  in  general 
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form  up  to  the  present  time,  and  recourse  must  be  had  in  most  cases  to 
qualitative  considerations,  approximate  solutions  and  computer  integration 
of  the  equations  of  motion .  The  research  reported  here  is  of  the  first  two 
types. 

C  )  RESULTS 

The  results  obtained  may  be  classified  according  to  whether  they 
deal  with  particles  which  are  reflected  by  the  magnetic  field,  or  which 
penetrate  the  field  to  tho  downstream  side  of  the  magnet.  Before  dealing 
with  these  results,  two  paragraphs  are  presented  containing  a  brief  review 
of  the  notation  and  the  basic  theory  utilized,  and  an  outline  of  the  in¬ 
jection  problem . 

(A)  Review.  Calculations  are  carried  out  in  the  standard  circular 

cylindrical  coordinate  system  (  p  ,  Q  ,  C)  =  (—- /  6  ,  — ),  where  a 

a  a 

is  the  field  extent  parameter  defined  above .  The  Lagrangian  for  an 
axially  symmetric  magnetostatic  system  maybe  written  as: 

A  =  y  (  b2  +  p2  &  +{2  )  +p©  «(p>,  C)  (2) 

whore  w(p,  £  )  is  the  dimensionless  magnetic  vector  potential  and 
time  measured  in  units  of  (  m/e  B0  ) .  The  system  is  conservative,  hence 
possesses  the  energy  integral 

p*  +  p2  e2  + 12  =  u2  ( 3 ) 
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(v  Is  Hie  dimensionless  velocity),  and  is  axially  symmetric  hence  pos¬ 
sesses  the  angular  momentum  integral 

p2e  +  p«(p,c)  =  Y  (4) 

The  latter  two  equations  imply  the  existence  of  forbidden  zones  whose 
boundaries  are  defined  by  the  equation 

y/p  -a(p  ,  C  )  *  +  u  (5) 

and  beyond  which  the  particles  of  given  (  u  ,  y)  cannot  penetrate. 

By  restricting  consideration  to  paraxial  motion  the  preceding  equation 
may  be  solved  for  p  to  give: 


p  +++  -  u  /  r0 

-1  *  ■Ji+SyVu*  '  ] 

(6) 

p  .++  =  u/ 

1  +  \J  i  +  2  y  rQ  /  -u2 

(7) 

p  ~+  =u/ro 

1  ♦  _ 

(8) 

P—  =  v/ro  1 

(9) 

where  ro  is  the  axial  variation  of  the  magnetic  field  magnitude. 

For  the  bell -shaped  field 

rQ  =  i/i  +  c2  do) 

(  b  )  Injection  Problem.  During  the  course  of  the  contract  it  became 
apparent  that  the  main  theoretical  problem  vould  be  to  predict  those 
values  of  the  appropriate  injoction  parameters  which  would  maximize  the 
radiation  emitted  in  a  given  spectral  region.  This  involves  defining  the 
proper  values  and  ranges  of  values  (allowable  sproads  )  of  the  coordinates, 
velocities,  principal  magnetic  field  value  (  Bq  )  and  inter- particle  phases, 
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with  the  emphasis  being  placed  on  determining  tho  valuos  of  the  coordinates 
and  velocities  as  the  most  important  parameters.  The  principal  results  of 
this  phase  of  the  investigation  are  to  be  found  in  the  a  fore- mentioned 
paper  and  will  be  quotod  briefly  below. 

(  c  )  Reflected  Particles.  For  those  particlos  which  do  not  penetrate 
to  tho  equatorial  plane  but  are  reflected  at  some  point  (p  ,  £  ), 
some  interesting  relations  may  be  derived.  The  following  expression  for 
the  principal  gyramagnotic  wavelength,  Xp  ,  maybe  derived  for 
paraxial  motion  in  a  relatively  strong  magnetic  field  assuming  that 
most  of  the  radiation  is  produced  near  a  reflection  point  where  the 
trajectory  may  be  approximated  well  by  a  circular  helix: 

xp  =  [°-0147  Bo  ro  ( Cr )  i-5TT  +  i  (11) 

where  V  is  the  particle  accelerating  voltage,  and  (  xp  >  B0  '  V  )  are 
moasured  in  (mm,  Kilogauss,  kV),  respectively.  Since  Xp  is  usually 
fixed  at  the  outset,  it  is  more  instructive  to  modify  this  equation  to 
obtain: 

V-511  0.0147  Bq  rQ  <Cr)X-l  (12) 

i.  K 

In  this  manner,  tho  concept  of  the  threshold  field,  B  ,  comes  into 

ot 

existence,  defined  to  be  tho  principal  magnetic  field  which  wi  II  just 
allow  radiation  of  wavelength  Xp  to  be  produced  at  the  reflection  point. 

This  relation  gives  the  first  formal  hint  that  tho  magnetic  field  magnitude 
is  the  main  limiting  systom  parameter  with  respect  to  high  radiation  output. 
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The  power  radiated  at  tho  reflection  point  may  be  shown  to  be 

P  =  Bo  ^o  (  «r>  [<  °-0147  )2  Bo  (  Cr)  "I  ]  (  13  ) 

where  P  is  measured  in  ev/ms  for  each  particle.  It  is  this  relation  which 
brings  out  most  strongly  the  importance  of  malting  the  magnetic  field  as 
high  as  possible;  the  gain  in  power  at  the  reflection  point  for  a  given 
wavelength  us  largo  as  BQ  goes  much  beyond  Bq^.  It  is  important  to  point 
out  also  that  the  voltage  which  must  match  Bq  through  the  second  previous 
equation  above  must  only  increase  in  a  linear  fashion  to  keep  pace  with  BQ, 
hence  tho  magnetic  field  is  again  shown  to  be  the  most  critical  system  pa¬ 
rameter.  There  exists  an  unknown  factor  in  this  analysis,  however.  It  is 
not  yet  known  whether  the  voltage  specified  by  the  equation  above  will 
be  sufficient  to  allow  tho  particle  to  penetrate  to  the  particular  reflection 
point,  assuming  optimum  values  for  the  other  injection  parameters.  In 
order  to  evaluate  this  factor  tho  reflection  point  parameters  must  be  re¬ 
lated  to  the  injection  point  parameters  through  the  solution  of  the  equations 
of  motion .  The  solution  must,  however,  be  a  higher  order  approximation  to 
the  trajectories  than  tho  paraxial  approximation,  bocause  no  reflection  is 
obtained  in  paraxial  (  first-order )  approximation.  A  reformulation  of  the 
equations  of  motion  has  been  carried  out  to  assist  this  effort.  We  introduce 
the  generalized  coordinates  (  i,  x)  the  meridional  plane,  where  if  is 
the  magnetic  scalar  potential  and  x  •*  the  magnetic  line  of  force  function. 

If  the  time  scale  is  implicitly  transformed  through  the  definition 

da/dT  =  T(p,  0  (14) 
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then  the  equations  of  motion  take  the  form: 

.  1  x*  SI'  1  do  .  *2 

*-  -  +  -  (u2- - )  - +  _  -  ( u2  -  if2  -  2  - )■  0  (15) 

r  r  p2  a*  P  at  p2 

v  r  i  ar  i  ap 

X  - -  +  —  (U2P2  -x2)  -  + - (x2  -  2x -x2) +  V  *  0  (16) 

r  r  ax  p  <sx 


where  the  dot  now  denotes  differentiation  with  respect  to  a  .  Although  a 
simple  first  order  solution  to  this  system  has  been  obtained,  no  solution 
has  yet  been  worked  out  which  allows  the  reflection  point  to  be  related 
to  the  injection  parameters.  This  constitutes  one  of  the  promising  directions 
for  further  research . 

(d)  Penetrating  Particles.  Those  particles  which  attain  the  equatorial 
plane  with  non-zero  forward  velocity  cannot  be  reflected,  hence  penetrate 
the  entire  field.  Their  zeroth  order  archetype  trajectory  is  the  paraxial 
trajectory,  v^ich  for  the  bell -shaped  field  has  the  form: 


o  =  c  i  1  +  r,5 


sin 


/i7 


1/4  u2  cot  Ca  + 


where  (  c,  ^> )  are  constants  to  be  evaluated  in  terms  of  the  injection 
parameters.  Based  upon  analysis  of  the  properties  of  the  forbidden 
zones  and  paraxial  trajectories  carried  out  in  the  a  fore -mentioned 
paper,  it  is  possible  to  make  certain  suggestions  concerning  desirable 
values  of  the  injection  parameters. 


(17) 
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D)  CONCLUSIONS 

In  order  for  the  particles  to  penetrate  to  the  high  field  regions  where 
the  magnetic  field  may  exert  a  sizeable  effect  it  seems  desirable  to  inject 
with:  ( a  )  negative  angular  momentum  constant  (  y  -  -  y  );(b)  zero 
angular  velocity  (9  =  0  )  ;  (  c  )  the  velocity  vector  initially  pointed 

along  a  line  of  force  {  o2  =  2y(  1  +  (}  ),  and  (  d  )  relatively  small  values 

o  ° 

mv 

of  )• 

Although  the  basic  elements  of  the  injection  problem  have  been  derived 
in  a  form  appropriate  to  guide  an  experimental  program,  a  great  deal  of  work 
remains  to  be  done.  More  exact  trajectories  will  have  to  be  obtained,  which 
eventually  will  have  to  include  the  radiation  reaction  effects.  Coherence 
conditions,  which  have  been  largely  ignored  during  the  present  contract, 
will  have  to  be  investigated.  These  problems  and  a  number  of  others  of 
varying  magnitudes  can  now  be  attacked  more  easily  and  intelligently  as 
a  result  of  the  work  reported  on  above  . 

Ill  )  EXPERIMENTAL  STUDIES 

The  experimental  program  carried  out  under  the  scope  of  this  project 
was  divided  into  two  distinct  phases: 

a  )  the  injection  studios  of  high-speed  electrons  through  strong  magnetic 
fields; 

b  )  the  studies  on  production  and  detection  of  cyclotron  electro¬ 
magnetic  radiation. 


•  10  - 


ADVANCED  KINETICS,  INC. 


A  )  EXPERIMENTAL  SET-UP: 

The  experimental  equipment  constructed  consisted  of  a  pulsed  electron 
source,  a  pulsed  high-field  magnet,  the  vacuum  system,  the  synchronizing 
electronics  and  the  detection  apparatus. 

i)  the  electron  sources  used  were  of  two  t/pes;  the  cold  cathode  type 
was  used  extensively  with  a  discharge  sustained  in  argon  by  application  of 
voltages  up  to  1  KV and  currents  in  excess  of  500  mA  .  Beams  of  several 
mA  were  extracted  up  to  10  KV  under  DC  conditions.  Alternatively,  the 
source  could  be  pulsed  by  triggering  a  high-power  line  discharge  circuit, 

in  which  case  beams  of  up  to  30  KV  and  about  10  mA  were  obtained  for 
times  of  the  order  of  2ps.  The  beams  obtained  have  a  diameter  of  about. 2  mm, 
on  impact  onto  a  quartz  plate  located  in  the  Faraday  cup,  after  traveling  through 
the  magnet  structure .  Hot-filament  cathodes  were  used  alternatively  to  obtain 
beams  of  equivalent  diameters  and  currents; 

ii)  the  pulsed  high  field  magnet  was  a  multi-turn  structure  similar  to  the 
onos  constructed  in  the  past  (  see  H  .  Furth  and  R  .  W.  Waniek,  Rev.  Sci .  Instr. 
27,  195,  (1956) .  Fields  up  to  250  kOe  could  be  obtained  over  1"  diameter  with 
a  field  extent  parameter  of  0.5"  .  After  the  injection  phase  was  completed,  this 
structure  was  modified  to  a  Helmholtz  geometry  of  equivalent  field  strength  in 
order  to  facilitate  the  insertion  of  hardware  close  to  the  high  field  region. 
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ill)  the  vacuum  system  was  constructed  so  as  to  allow  under  usual  in¬ 
jection  conditions  a  vacuum  of  around  I0“®  Torr.  This  pressure  could  be 
increased  to  permit  a  visible  observation  of  beams  injected  under  D.C  . 
conditions,  The  vacuum  in  the  cold-cathode  electron  source  was  around 
100  microns  and  was  regulated  by  a  feed  valve.  All  probes  were  fed 
through  vacuum-tight  seals. 

iv)  the  synchronizing  electronics  allowed  a  precise  timing  of  the  in¬ 
jection  pulse  with  respect  to  the  high  field  magnet  operation.  Under 
transient  conditions  the  electrons  were  injected  close  to  the  maximum 
of  the  field  pulse  to  avoid  extraneous  fi  or  LI  effects  in  the  system. 

This  was  about  70  ps  after  zero  time,  the  magnet  pulse  lasting  close 
to  140  ps.  A  master  signal  was  used  to  trigger  the  magnet  circuitry, 
and  after  the  selected  time-delay,  the  high  voltage  line  was  discharged 
across  the  accelerating  electrode. 

v)  Two  detection  systems  were  adopted,  one  for  each  phase  of  the 
experiment: 

a)  during  the  injection  phase  a  specially  shielded  Faraday  cup 
was  developed  with  triaxial  vacuum-tight  connection.  The  beam  current 
was  monitored  on  a  CRO  by  displaying  the  voltage  pick-up  across  a 
4.7  KQ  resistor  joining  the  cup  to  ground.  The  cup  could  be  moved 
axially  about  8",  allowing  to  explore  a  region  extending  from  5"  below 
the  magnet  to  the  top  of  the  magnet.  All  CRO  readings  were  made 
inside  a  special  screened  cage. 
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b)  the  radiation  detection  phase  included  first  an  attempt  to  monitor 


possible  signals  with  220  kMe  horn  and  crystal  unit.  There  was  considerable 
doubt  that  this  unit  could  detect  the  radiation  because  of  unknown  sensi¬ 
tivity  of  the  crystal  in  the  region  beyond  300  kMc  (  =  1  mm  ) ,  and  because 
of  the  unknown  response  time.  Experiments  were  conducted  to  test,  and 
possibly  to  improve,  sensitivity  and  response  times.  IN23  and  MA  -408  -B 
crystals  were  subjected  to  pulsed  signals  from  an  X-band  magnetron  system 
(  2J51  )  operating  around  19  kMc  .  Different  conditions  of  reverse  bias  were 
attempted.  It  was  noticed  that  crystals  with  good  sensitivity  had  good  rise¬ 
time,  and  vice-versa  .  Crystals  with  good  rise- time  and  sensitivity,  however, 
were  generally  unaffected  by  the  biasing  procedure,  whereas  poor  crystals 
improved  such  properties  upon  biasing.  These  results  together  with  the 
consideration  of  the  fragility  of  the  crystals  used  at  G-band  indicated  that 
biasing  was  not  advisable  in  this  operation.  Rise- times  of  less  than  50 
nanoseconds  have  been  measured  with  a  Fischer  source  directed  at  the 
horn . 

Since  the  power  expected  from  the  radiator  decreases  at  longer 
wavelength  (  S  =  v2  B2  and  w  =  eB/my  )  a  critical  problem  in  this  experiment 
was  the  development  of  a  detection  system  operating  with  good  efficiency 
in  the  region  below  500  microns  for  pulsed  signals. 

During  the  last  quarter  a  high  sensitivity  PEM  detector  was  developed 
to  be  operated  in  liquid  helium  at  fields  around  10  kOe.  Such  system  was 
tested  in  helium  with  superconducting  fields  up  to  24  kOe  but  has  not  yet 
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been  connected  to  the  radiator  itself.  This  system  is  expected  to  be 
sensitive  in  the  region  of  interest. 

B)  RESULTS 

a)  INJECTION  EXPERIMENTS: 

The  transmission  of  electrons  through  the  strong  magnetic  fields  was 
studied  as  a  function  of  injection  energy  and  angle .  The  strength  of  the 
magnetic  field  was  varied  stepwise  to  fields  beyond  100,000  oersteds. 

The  injection  was  performed  on-axis  and  off-axis.  Also  various  positions 
of  the  detecting  Faraday  cup  and  of  the  injecting  electron  source  were 
used  .  The  data  had  to  be  corrected  to  eliminate  the  effect  of  stray 
electrons  (  specular  reflection  )  by  evaluating  the  Larmor  radii  of  electron 
orbits  belonging  to  each  particular  injection  angle.  When  such  Larmor  radii 
were  smaller  than  the  wall  dimension,  the  particular  orbit  was  accepted  as 
falling  within  the  limits  of  confidence  .  Orbits  touching  the  walls  before 
passage  through  the  magnot  were  rejected. 

The  effect  of  the  forbidden  zones  on  the  transmission  of  the  electron 
beams  through  the  strong  field  was  noticeable.  For  instance,  for  20  KeV 
injection,  the  intensity  between  0°  and  3.5°  injection  would  vary  by  more 
than  a  factor  of  2 .5  (  Fig.  I).  Similar  conditions  are  obta ined  also  for  the 
off-axis  injection  (  Fig.  2  ).  The  strong  magnetic  fields  increase  considerably 
the  transmission .  At  fields  around  20  kOe  the  current  intensity  drops  by  more 
than  a  factor  of  three,  at  a  constant  angle,  with  respect  to  the  zero  angle 
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injection, and  fields  in  excess  of  100  kOe  bring  back  the  intensity  to  the 

original  0°  injection  value.  All  the  above  injection  experiments  were 

carried  out  with  nearly  zero  (  to  within  the  setting  accuracy  of  the 

electron  source  )  angular  momentum. 

The  results  show  definitely,  therefore,  that  the  forbidden  zones  impair 

electron  injection  ( roughly  a  factor  of  5  at  0°  and  100  kOe  )  but  that 

transmission  channels  existed  for  the  geometry  studied  at  around  3.5° . 

At  such  injection  angle  and  under  the  application  of  fields  around  100,000 

Oe,  currents  of  up  to  30  pA  could  be  transmitted  corresponding  thus  to  a 
0 

2  ps  burst  of  3 .6 '10  electrons.  For  a  10  KeV  electron  bunch  at  100  kOe 

-14 

the  cyclotron  radiation  at  roughly  300  kMc  would  be  about  8 '10  watts/ 
electron  yielding  a  power  of  about  6  pW  into  the  solid  angle  used. 

B  )  DETECTION  STUDIES 

With  an  estimated  crystal  sensitivity  of  about  1  mV/pW  at  300  kMc, 
the  experiment  with  the  present  electron  currents  is  just  marginal  and  the 
outcome  would  clearly  depend  on  problems  of  amplification .  The  220  kMc 
horn  and  detector  were  positioned  in  a  meridional  plane  outside  the  magnet 
and  aligned  with  the  center  of  the  same .  Special  precautions  had  to  be 
taken  to  avoid  the  strong  level  of  noise  pick-up  during  the  pulse  phase. 

The  signal  from  the  detector  was  fed  via  a  175  Ohm  cable  to  a  system  of 
2  cascaded  Tektronix  preamplifiers,  type  L.  The  maximum  sensitivity  of 
such  a  system  operated  with  appropriate  termination  for  optimum  noise 
reduction  was  1  mV/cm. 
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Under  these  conditions  results  were  obtained  during  one  run  indicating 
the  following  signals  for  injection  at  5°: 


Magnetic  Field 

Peak  Signal 

Cyclotron  Frequency 

0  KOe 

20  mV 

(  background  due  to  induced 

effects  of  injected  electron  pul 

87  KOe 

28  mV 

250  kMc 

108  KOe 

42  r.  '.V 

310  kMc 

Within  the  experimental  error  (  estimated  at  about  10%  ),  the  signal  would 

seem  to  indicate  a  square  power- field  relationship,  as  expected  for  cyclotron 

radiation  production.  The  signal  obtained  at  108  kOe  would  seem  to  be  clearly 

above  the  noise  level  by  a  factor  of  two  (  at  85  kOe,  the  expected  power  is 
-14  -13 

4*10  w/electron  versus  about  10  w/electron  at  108  kOe  ) .  Under  the 
above  assumption  of  a  crystal  sensitivity  of  1  mV/pW,  such  signal  would  corre¬ 
spond  to  some  20  pW  of  radiated  powor  which  is  about  a  factor  of  3  higher  than 
the  theoretical  figure.  Such  factor,  however,  might  well  be  attributed  to  the 
crystal  sensitivity. 

CONCLUSION 

The  evidence  seems  to  indicate  that  20  pW  of  power  have  been  extracted 
at  300kMc  from  a  bunch  of  10  keV  electrons  by  cyclotron  radiation  in  a  field 
of  108,000  Oersteds.  The  presently  unknown  sensitivity  of  the  detector  used 
could  reduce  the  power  figure  to  the  theoretically  expected  of  6  pW .  The 
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encouraging  evidence  of  the  above  experiments  indicates  that  transmission 
and  radiation  can  be  achieved  in  high  field  geometries.  The  extension  of 
this  work  must  clearly  Involve  higher  injection  velocities,  higher  fields 
and  the  necessary  radiation  detector  for  the  shorter  wavelengths.  The 
injection  at  20  keV  through  200,000  Oe  is  expected  to  increase  the  power 
by  a  factor  of  1b  to  some  200  pW  radiated  at  550  microns  or  545  kMc . 
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"THE  INJECTION  OF  HIGH  VELOCITY  CHARGED  PARTICLES 


INTO  STRONG  MAGNETIC  FIELDS"* 
by 

Malcolm  A .  MacLeod  and  Ralph  W.  Waniek 
Advanced  Kinetics,  Inc. 

Costa  Mesa,  California 

ABSTRACT 

Strong  magnetic  fields  are  useful  in  work  with  high  velocity  charged  particles 
by  virtue  of  their  ability  to  rotate  the  particle  momentum  vector.  The  degree  of 
usefulness  is  measured  by  the  extent  to  which  the  desired  rotations  can  be  produced 
at  the  proper  positions  in  the  field.  It  is  frequently  necessary  to  inject  charged 
particles  into  magnetic  fields  at  points  remote  from  these  proper  positions,  in  which 
case  the  injection  parameters  which  will  enable  the  particles  to  reach  these  posi¬ 
tions  must  be  specified .  For  the  production  of  gyromagnetic  radiation  in  axially 
symmetric  magnetostatic  fields  the  proper  positions  correspond  to  the  high  field 
regions  near  the  symmetry  axis.  The  present  introductory  study  is  concerned  with 
delineating  approximate  injection  parameter  ranges  which  will  allow  the  parti¬ 
cles  to  attain  the  regions  of  high  magnetic  field  strength  in  a  field  whose  axial 
variation  is  bell-shaped.  The  single-particle  model  is  adopted,  and  the  results 
of  both  qualitative  and  quantitative  investigations  are  presented. 

•Research  supported  by  the  Electronics  Research  Directorate  of  the  Air  Force 
Cambridge  Research  Laboratories. 
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A  general  formulation  of  the  axially  symmetric  magnetostatic  field  forbidden 
zone  theory  is  given  and  the  useful  concept  of  the  paraxial  forbidden  zones  is 
introduced.  This  theory  is  then  applied  to  obtain  qualitative  answers  to  the  in¬ 
jection  problem  for  the  bell-shaped  field  for  both  zero  and  non-zero  values  of 
the  angular  momentum  constant.  Quantitative  relations  for  the  zero  angular 
momentum  constant  case  are  obtained  using  the  bell-shaped  field  paraxial 
trajectories.  A  discussion  of  the  limitations  of  the  approximations  employed 
is  included . 
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I  .)  INTRODUCTION 

The  recent  development  of  strong  pulsed  magnetic  fields  in  the  megagauss 
region  '  hai  made  possible  a  number  of  interesting  applications,  such  as  co¬ 
herent  electromagnetic  radiation  production  and  momentum  selection  of  high 
energy  charged  particles  .  Many  of  these  applications  involve  the  introduction 
of  a  beam  or  bunch  of  particles  into  a  pre-existing  magnetic  field.  One  of  the 
central  problems  in  such  work  is  the  choice  of  beam  injection  parameters  which 
will  optimize  the  desired  behavior  of  the  particles  in  the  field.  This  paper  is 
concerned  with  the  zeroth  order  delineation  of  those  ranges  of  injection  pa¬ 
rameters  for  a  particular  type  of  field  configuration  which  will  allow  the  parti¬ 
cles  to  attain  regions  of  high  magnetic  field  strength . 

The  problem  in  the  forefront  during  the  course  of  this  work  has  been  the 
production  of  gyromagnetic  radiation,  i  .e  .,  the  production  of  electromagnetic 
radiation  by  the  radial  acceleration  of  charged  particles  in  a  magnetic  field. 
This  is  essentially  a  relativistic  effect  whereas  the  mathematical  model  of  the 
physical  situation  presented  here  is  of  the  nonrelativistic  variety,  hence,  due 
caution  must  be  exercised  in  tho  interpretation  of  results  and  their  application 
to  the  study  of  gyromagnetic  radiation  production . 

In  the  interest  of  broad  applicability  of  the  results,  the  current  paper  deals 
only  with  the  choice  of  the  particle  injection  parameters  necessary  to  attain 
high  field  regions  in  a  given  magnetic  field  configuration.  The  choice  of  the 
field  parameters  is  determined  by  the  particular  application  and  therefore  falls 
outside  the  scope  of  the  present  discussion. 
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II .)  GENERAL  CONSIDERATIONS 
A.)  THE  MODEL 

Two  aspects  of  the  modeling  problem  exist:  the  choice  of  a  representative 
physical  system  and  the  degree  to  which  the  mathematics  of  this  system  must  be 
simplified  to  obtain  tractable  equations  which  yield  meaningful  solutions. 

Thus  the  problem  arises  of  the  choice  of  a  representative  physical  system 
which  is  complicated  enough  to  enable  non-trivial  results  to  be  extracted  while 
not  being  too  complicated  so  as  to  impede  the  analytical  development.  Fortu¬ 
nately,  such  a  system  is  suggested  by  strong-field  magnet  configurations. 

In  the  first  place,  t;nical  strong-field  magnets  possess,  at  least  approxi¬ 
mately,  an  axis  of  symmetry,  hence  the  magnetic  field  may  be  taken  to  be 
axially  symmetric . 

Secondly,  for  pulsed  fields  the  ratio  of  the  period  of  the  field  to  the 
period  of  time  spent  by  the  particles  in  the  field  may  frequently  be  assumed 

3 

to  be  large.  Typical  field  periods  are  1  -  10  psec,  whereas  the  longest 
particle  periods  are:  0.05  psec  (  1  Kev  electron  )  and  2.0  psec  (  1  Kev  proton  ), 
both  particles  having  a  Ira  path  length  in  the  field.  Thus  the  field  may 
usually  be  taken  to  be  magnetostatic  to  a  high  degree  of  approximation . 

Thirdly,  pulsed  high  field  magnets  concentrate  their  fields  within  certain 
regions  along  the  symmetry  axis,  frequently  displaying  a  symmetry  about  an 
"equatorial"  plane  perpendicular  to  the  axis  at  the  magnet  center.  In  fact, 
it  is  the  spatial  inhomegeneity  thus  produced  that  creates  the  injection  prob¬ 
lem.  In  such  configurations,  the  field  magnitude  is  zero  at  +  “along  the 
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axis  and  gradually  rises  to  a  local  maximum  at  the  magnet  center .  This  behavior 
is  adequately  represented  for  the  present  purposes  by  a  magnetic  field  of  the  type 
whose  axial  dependence  is  given  by  the  "bell -shaped"  curve. 

B  (  0,  0,  z  )  =  B0  /  1  +  (z/a)2  (  1  ) 

*  axis  is  the  axis  of  symmetry  and  the  a  denotes  the  half-width 
of  the  field  intensity  rui ve  at  I  kj  If- maximum .  This  field  has  the  advantage  that 
many  of  its  properties  are  amenable  to  analytic  treatment^. 

Lastly,  the  ratio  of  the  mean  free  path  of  the  particles  in  the  beam  or  bunch 
to  the  local  gyromagnetic  radius  will  be  assumed  to  be  large;  space  charge  effects 
will  be  ignored.  Roughly  speaking,  for  instantaneous  particle  currents  of  a  few 
milliamperes  and  kinetic  energies  of  at  least  a  few  Kev  there  is  no  space  charge 
problem.  To  verify  this  remark  and  to  estimate  these  effects  more  closely,  con¬ 
sult  Lopukhin^.  Here  the  beam  will  be  taken  to  be  an  aggregate  of  single  par¬ 
ticles. 

The  simplifications  necessary  to  establish  a  tractable  mathematical  model 
are  considerably  more  restrictive  in  nature  than  those  introduced  above.  Since 
the  exact  differential  equations  of  motion  of  a  classical  structureless  point 
charged  particle  moving  in  a  magnetostatic  field  have  never  been  integrated 
completely^,  one  would  expect  to  introduce  some  simplification  to  obtain 
analytic  results.  The  main  assumptions  introduced  here  are  two:  radiation 
reaction  forces  may  be  neglected,  and  only  the  paraxial  motion  will  be  con¬ 
sidered  . 
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The  neglect  of  the  radiation  reaction  forces  actually  involves  two  «epa ra te  con¬ 
ditions,  since  there  are  both  relativistic  and  nonrelativistic  contributions  to  these 
forces.  By  restricting  consideration  to  nonrelativistic  motion  (  v/c  *  8  <<1  )  the 
exact  differential  equations  of  motion  reduce  to  the  following^: 

dv/dt  -  e/m  v  xB  =^l/d)^d2v/dt2^  (2) 

where 

1/d  =  e2/6*fi0mc^  (3) 

and  where  the  right-hand  side  is  the  well-known  nonrelativistic  contribution  to 
the  radiation  force.  Further  analysis  shows  that  this  contribution  may  be  neglected 
when^ 

a2  «  dAp  (4) 

where  <u  is  the  gyromagnetic  (  cyclotron  )  frequency,  Tp  is  the  time  which  the  par¬ 
ticle  spends  in  the  field.  By  using  the  definition  of  the  gyrofrequency 

jo  =  eB/m  (  5  ) 

the  inequality  (  4  )  may  also  be  written  as 

B2  *"<  (  m/e  )2  dA  (  6  ) 

P 

This  restriction  may  be  shown  to  be  mild  for  the  type  of  motion  considered  here, 
though  the  restriction  is  more  severe  for  particles  of  low  mass  than  for  particles 
of  high  mass.  It  should  be  noted  that  longer  path  lengths,  such  as  would  be  en¬ 
countered  in  a  mirror  device,  could  make  the  neglect  of  the  radiation  reaction 
force  become  a  serious  error. 
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The  second  simplifying  assumption  leads  to  the  analysis  of  paraxial  motion. 

The  restrictions  involved  in  this  assumption  are  three: 

1)  the  particle  trajectory  radius  ( the  deviation  of  the  trajectory  from  the 
field  axis )  must  be  small  compared  to  a  characteristic  dimension  of  the 
field  source,  denoted  by  a 

r  «  a  (7) 

2)  the  slope  of  the  trajectory  must  be  small: 

dr/dz  «  1  (8) 

and 

3)  the  ratio  of  the  gyromagnetic  radius  to  a  characteristic  field  length 
must  be  large  compared  to  the  vector  potential 

mv/eaB  >>  Ag  (  9 ) 

The  effect  of  these  restrictions  on  the  particle  motion  is  best  discussed  after 
the  trajectories  have  been  obtained. 

It  is  clear  that  the  particles  may  penetrate  to  high  field  regions  in  an  axially 
symmetric  magnetostatic  field  when  injected  on  the  axis  with  zero  slope  because 
the  Lorentz  force  is  zero  there  and  the  particles  pass  through  the  field  undeviated. 
Such  a  path  is  uninteresting  because  the  field  does  not  influence  the  particle  motion 
at  all .  A  particle  injected  off-axis  with  zero  slope  or  on-axis  with  a  small  slope 
will,  however,  be  acted  upon  by  the  magnetic  field  and  will  pass  through  the  field 
provided  certain  conditions  are  fulfilled.  It  is  the  approximate  delineation  of 
these  conditions  and  of  the  high  field  behavior  of  the  particles  that  is  the  object 
of  the  present  study,  and  it  is  here  that  the  paraxial  approximation  proves  useful . 
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The  physical-mathematical  model  is  now  complete.  The  results  presented 
below  are  based  on  the  study  of  the  nonrelativistic  paraxial  motion,  neglecting 
radiation  reaction,  of  a  single  charged  particle  moving  in  an  axially  symmetric, 
magnetostatic  field  whose  axial  variation  is  bell-shaped . 

B)  UNITS 

It  is  necessary  in  order  to  obtain  the  greatest  clarity  of  exposition  to  adopt 
a  simple,  yet  physical,  system  of  units  for  the  variables  involved  in  the  theory. 

The  system  studied  requires  for  its  complete  characterization  the  specification 
of  eleven  parameters  which  may  be  taken  as  follows:  the  charge  e  and  the 
mass  m  of  the  particles,  tho  injection  point  values  of  the  variables  (  q^,  cjj0, 
t  (  i  -  1,  2,  3  ) ),  a  characteristic  field  value  and  a  field  extent  parameter. 
Although  a  more  complex  fiold  may  involve  any  numbor  of  spatial  parameters, 
the  bell-shaped  field  involves  only  one,  which  is  taken  to  be  a,  the  half-width 
at  half-maximum  of  the  axial  fiuldcurve.  The  logical  choice  for  the  character¬ 
istic  field  value  is  then  B  .  tho  maximum  axial  field  value. 

o 

The  dynamical  formulation  may  be  completely  freed  of  its  explicit  dependence 
on  the  system  parameters  by  an  appropriate  choice  of  units.  The  choice  made  here 
is  as  follows: 

Length:  a 

Magnetic  Field:  BQ 

Time:  b -■  m/eB0  (10) 
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Note  that  time  Is  now  measured  in  units  of  the  reciprocal  gyromagnetic 
frequency  with  respect  to  the  field  center. 

The  natural  coordinate  system  to  employ  for  the  investigation  of  a 
general  axially  symmetric  field  is  the  standard  circular  cylindrical  system 
( r,  6  ,  z  ),  or  (  p  ,  6  ,  C  )  in  terms  of  the  dimensionless  coordinates  defined 
by  the  equations*: 

ps  r/a 

'  =■  z  /a  (11) 

*  Dimensionless  quantities  will  consistently  bo  denoted  by  Greek  Symbols. 

Thus  the  dimensionless  time  and  velocity  take  the  form 

t  =  t/b 

u  =  m  v  /  eaB0  (  12  ) 

while  the  magnetic  field  and  vector  potential  become: 

r  =  B/^0 

a  =  A/aB0  (  13) 

A  very  convenient  aspect  of  the  above  definitions  is  that  the  magnitude 
of  the  dimensionless  velocity,  u  ,  is  in  fact  the  dimensionless  gyro- 
magnetic  radius  of  the  particle  at  the  field  center,  assuming  the  velocity 
vector  to  be  entirely  in  the  equatorial  plane. 

u  -  mv/oaBQ  (  14  ) 

This  quantity  is  one  of  the  two  most  important  trajectory  parameters.  It  is  here 

that  the  advantage  of  strong  magnetic  fields  is  apparent,  for  high  velocities  may 

be  used  and  u  may  be  still  kept  relatively  small,  allowing  the  influence  of  the 
magnetic  field  to  be.  felt. 
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III.)  FORMULATION  OF  THEORY 
A.)  FORBIDDEN  ZONES 

The  concept  of  the  forbidden  zone*  was  introduced  by  C .  Stdrmer  (  7  )  in 
connection  with  the  problem  of  the  motion  of  a  charged  particle  in  the  geo¬ 
magnetic  field.  The  following  generalized  treatment  present*  a  formulation 
tailored  to  the  problem  to  be  investigated. 

The  Lagrangian  point  of  view  will  be  adopted  here  as  the  basis  for  the 
dynamical  formulation.  The  nonrelativistic  Lagrangian  fora  single  charged 
particle  moving  in  a  magnetostatic  field  is  known  to  be: 

L  =  l/2mv*  +  ev  '  A  (  15  ) 

On  introducing  the  variables  defined  above  and  denoting  the  dimensionless 
Lagrangian  by: 

A  E  m  L/  (eaB^2  (16) 

the  expression  (  15  )  for  the  Lagrangian  bocomes: 

A  =  1/2  u2  +  u  •  a  (  17  ) 

Since  it  can  be  shown  that  the  magnetic  vector  potential  of  an  axially  symmetric 
magnetostatic  field  may  be  specified  completely  by  a  vector  with  only  an  azi¬ 
muthal  component,  the  basic  Lagrangian  takes  the  form: 

A  =  1/2  (p2  +  p2  02  +  £)  +  p  9  a  (  p  ,  '  )  (  18  ) 

It  is  evident  that  9  is  a  cyclic  variable,  hence  the  corresponding  momentum  is 
conserved: 

P 2  0+  P  or  (  p  ,  C )  =  y  =  const .  (  19  ) 
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Furthermore,  the  kinetic  energy  is  conserved  because  the  field  is  magneto¬ 
static: 


p2  +  p2  d2  +  £*  *  U2  =  const. 

(  Note  that  the  kinetic  energy  is  defined  here  to  be  one-half  of  the  left- 
hand  side  of  equation  (  20  ) ) .  Eliminating  0  between  the  preceding  two 
equations,  the  following  expression  is  obtained: 

fJ2  +  C2  +  (“  -  a)2  =  u2 

P 

Clearly 

1  v2  J 

(  P  -  «)z  <  v/ 

Note  that  when  the  equality  sign  holds  in  (  22  )  the  p  and  C  velocities 
are  zero,  the  velocity  is  entirely  azimuthal,  and  a  two-parameter  family 
of  curves  is  defined  in  the  meridional  plane: 

y/p  -  or  (  p  ,  C  )  *  t 

These  curves  are  the  boundaries  of  the  forbidden  zones;  they  delineate 
those  regions  of  the  (  o,  f  )  plane  which  are  inaccessible  to  particles 
characterized  by  the  parameters  (  u,  y  ) . 


(20) 


(21) 

(22) 


(23) 


B)  PARAXIAL  FORBIDDEN  ZONES 

Although  equation  (  23  )  defines  the  forbidden  zones,  it  does  so 

implicitly  for  any  field  once  the  axial  field  variation  is  known  .  In  spite 

of  the  fact  that  the  vector  potential  for  the  bell-shaped  field  is  known  in 

4 

closed  form  (  see  Glaser  )  and  that  the  exact  zones  may  be  obtained 

8 

analytically  ,  it  is  convenient  to  limit  the  present  discussion  to  a  con¬ 
sideration  of  the  paraxial  zones. 
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It  is  well  known  that  the  vector  potential  of  an  axially  symmetric 

magnetostatic  field  may  be  represented  as  a  series  expansion  about  the 
.4 

symmetry  ax  is  :  05 

p  £  (“I)"  flllSl  P  2n 

a?  (p,  r)  =  "j"  n  =  0  n'.  (n  +  1)'.  2n  (7  )  ( 24 ) 

=  "2  l^r(  c)  -P2/8  (  d2  r(  r)/d  r2)  +  ••••  j  (25) 
where  r (  C,  )  denotes  the  axial  variation  of  the  dimensionless  field  magnitude: 

HO  =  IA>  'B  *B  (26) 

If  the  axial  field  variation  is  reasonably  smooth  and 

P2  <<  1  (27) 

the  first  term  in  the  expansion  will  give  a  good  approximation  to  the 

vector  potential .  This  first  term  is  the  paraxial  vector  potential,  (p,  r): 

P 

or0  (p,  o  -  T  r(0  (28) 

and  the  paraxial  forbidden  zones  are  derived  by  substituting  expression  (  28  ) 
into  equation  (  23  ): 

v/p  ~\p/2  )r  (  ")  =  +  u  (29) 


The  explicit  form  of  the  zones  p  (  C.)  is  given  now  by  the  solutions  of  the 
quadratic  equation 

2 

r(')  p  +  ut-y=o  (30) 

j- 

Since  y  may  be  positive  or  negative  and  the  solution  to  the  quadratic 
introduces  two  more  sign  possibilities  besides  the  two  already  evident  in 
equation  (  30  ),  there  are  eight  roots,  which  may  be  denoted  by  oa|jCj 
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where  a  refer*  to  the  sign  of  the  first  degree  term  In  equation  (  30  ),  b 
refere  to  the  sign  of  y  and  c  refers  to  the  sign  introduced  by  the  form  of 
the  solution  to  the  quadratic.  After  excluding  the  unphysical  roots 
(  corresponding  to  p  negative  real  or  complex  ),  the  following  four 
roots  are  left: 


«  ! 

-  1  J]  +  2y  r  (  r  ) / u2  | 

(31  ) 

P+++  =  r(c)  j 

u 

r  (C) 

f  - — — . 

1  +  V  1  +  2  y  r  (  r  )/J 

1  - -1 

(32) 

; 

o  -«+  =  u  i 

r(c) 

V 

l  +  V  l  -  2  y  r  ( r  )/u2 

I  ( 33  ) 

p—  -  r(  c) 

1  -  v  1  -  2  V  r  ( r )/  u2 

1  (34) 

where  y  indicates  the  magnitude  of  y  .  Note  that  when  y  -  o,  p+++ 
and  p—  both  vanish,  whereas  o«++  and  p~+  coalesce  to  give  a  tingle 
zone  curve 

P0  =  2u/r(;)  (35) 

and  the  zones  are  seen  to  be  contour  lines  of  the  vector  potential .  Certain 
general  features  of  those  zones  will  now  be  pointed  out  briefly. 

It  is  clear  from  equations  (  31  )  and  (  32  )  that  solutions  are  possible 
for  all  values  of  £  (  for  all  values  of  f  (  £)  )  when  y  is  positive,  whereas 
solutions  for  negative  y  are  possible  only  for  those  values  of  r  for  which 

r(C)  <  u2/2y  (36) 

or  for 

r}  >  2y/  u2  -  1  ( 37 ) 
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If 


2y/u*  <  1 

then  both  branches  of  the  negative  y  forbidden  zone  exist  for  a  1 1  values 
of  f  because  the  inequality  (  37  )  is  trivially  true.  However,  if 

2v/u2  >  1 

then  there  exists  a  lower  limit  on  £  ,  given  by 


(38) 


(39) 


/ 

r  =  +  \l  2y/v 2  -  1  ( 40 ) 

for  which  the  forbidden  zone  exists.  At  this  point  the  equality  sign  holds 
in  (  36  )  and  < 


o — *■  =  p—  ( 41  ) 

Thus  it  becomes  clear  that  p  and  P___  are  simply  different  branches 

of  the  same  curve,  p _ +  always  being  larger  than  p _ ,  except  at  their 

point  of  intersection .  No  such  problem  arise  for  posititve  y  .  p(  ( ,  and 
P-++  exist  for  all  values  of  '  ,  they  never  intersect,  and 

o-++  >  o+++  (42) 


In  fact, 

o_++  -  P+++  =  2u/  r(0  (43) 

Now  as  £  -*  00 , 

2y  r  (  ; )/ u2  -o  (44) 

and  the  square  roots  in  equations  (  31  )  -  (  34  )  may  be  expanded  by  the 
binomial  theorem  to  give 

P+++  •*»*  y/w  (45) 
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p-++  •  2Vr<;) 

+  y/u 

(46) 

p~f  *2u/r(C) 

-  v/’J 

(47) 

P —  »  y/u 

(48) 

These  relations  show  that  for  large  £  the  roots  p_(_)_+  and  p__  are  governed 

largely  by  the  angular  momentum,  whereas  p_t  (  and  P _ f  are  governed 

largely  by  the  magnetic  field.  Since  the  paraxial  approximation  has  limits 
for  some  value  of  p  for  every  value  of  C  within  which  it  can  reasonably  be 
applied,  and  since  the  roots  values  given  by  equations  (  45  )  and  (  46  ) 
blow  up  as  ^  ~  “  ,  the  corresponding  zones  should  be  employed  with  discretion 
(  see  discussion  below  on  the  inaccuracies  introduced  by  employing  the  paraxial 
approximation) . 

Note  that  the  zones  given  by  equations  (  45  )  and  (  48  )  (  and  the 
corresponding  more  accurate  expressions  (  31  )  and  (  33  )  appear  only  as  a  result 
of  the  particle  having  non -zero  angular  momentum  about  the  symmetry  axis. 
These  zones  are  the  centrifugal  barrier  zones  which  arise  from  the  contribution 
of  the  centrifugal  potential  to  the  meridional  plane  potential  and  prohibit  the 
particle  from  reaching  the  axis.  The  other  zones,  given  by  equations  (  32  ) 
and  (  34  ),  arise  mainly  from  the  vector  potential  contribution  to  the  meridi¬ 
onal  plane  potential  which  acts  to  prevent  the  particle  from  approaching 
more  closely  to  the  field  source  than  its  kinetic  energy  will  allow.  (  The  closer 
the  field  source  is  approached,  the  higher  does  the  magnetic  field  become;  this 
action  forces  the  particle  away  from  the  high  field  regions  and  creates  the  major 
part  of  the  injection  problem) . 
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More  details  on  the  forbidden  zones  will  be  given  in  the  section  below 
which  deal  with  the  injection  conditions.  Attention  will  now  be  turned  to  the 
background  for  the  derivation  of  the  paraxial  trajectories. 


C  .)  THE  MERIDIONAL  PLANE  POTENTIAL 


It  is  well  known  that,  under  the  approximations  enunciated  above, 
the  motion  of  a  charged  particle  in  an  axially  symmetric  magnetostatic  field 
may  be  considered  to  take  place  in  the  meridional  plane  which  is  defined  by  the 
axis  of  symmetry  and  the  instantaneous  particle  position.  As  the  particle 
traverses  the  field  the  meridional  plane  rotates  with  the  particle  angular 
velocity.  The  motion  in  this  plane  moy  be  shown  to  be  conservative  and  to 
take  place  under  the  influence  of  a  meridional  plane  potential,  »,  which  is 
here  defined  as 


-2 


it  (  o ,  r  )  =  |/2 


y/o  -a  (  P  /  r) 


J 


(49) 


(  see  Glaser  ) .  This  decoupling  of  the  meridional  and  azimuthal  motion 
is  made  possible  by  the  existence  of  the  second  integral  of  the  motion,  the 
angular  momentum  integral  given  in  equation  (19)  above.  In  fact,  con¬ 
sidering  the  motion  as  taking  place  under  the  potential  n  is  just  a  question 
of  interpreting  equation  (  21)  above  in  the  proper  manner.  A  good  deal  of 
insight  into  the  types  of  particle  motion  may  be  gained  by  a  study  of  the 
meridional  plane  energy  relations. 

As  far  as  the  energy  representation  is  concerned,  the  particle  motion 
is  determined  by  the  parametors  (  u  ,  y)  and  their  associated  surfaces  in  the 
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(  p,  Q  ,  e  )  space  ( the  energy,  e,  is  visualized  as  being  plotted  vertically 
over  the  meridional  plane  ).  The  surface  associated  with  the  particle  velocity 
constant,  u  ,  is  a  plane  parallel  to  the  meridional  plane  at  e  *  u2/2  ,  while 
the  surface  associated  with  the  particle  angular  momentum  aonstont,  y  ,  is  the 
meridional  plane  potential  energy  surface  e  =  it  (  p  ,  Q  ).  The  magnitude 
of  the  meridional  plane  kinetic  energy  is  determined  uniquely  for  each  point 
(  P  /  Q  )  by  the  relation 


1/2  (p2  +  C2)  =  1/2  j^u2-(  t  -a(p, 
hence,  it  is  simply  equal  to  the  vertical  distance  between  the  two  surfaces. 

Since  the  forbidden  zone  boundaries  are  defined  to  be  the  aggregate  of  those 
points  for  which  the  meridional  plane  kinetic  energy  is  zero,  they  are  now  seen 
from  a  geometrical  point  of  view  as  the  curves  which  constitute  the  intersection 
of  the  u  plane  with  the  tr  surface . 

There  are  two  mechanisms  represented  in  the  form  of  * .  One  is  that 
associated  with  the  magnetic  field  which  gives  rise  to  the  presence  of  the 
vector  potential  in  it.  The  other  is  that  associated  with  the  conservation  of 
angular  momentum  which  gives  rise  to  the  centrifugal  potential  contribution 
to  tr.  Asa  particle  tries  to  approach  either  the  field  source  or  the  symmetry 
axis,  it  encounters  strong  forces  which  turn  it  aside,  as  shown  by  the  increasing 
slope  of  the  potential .  The  distance  of  closest  approach  to  either  region  depends 
mainly  on  the  value  of  the  velocity  constant,  u,  for  a  given  field  and  y  value. 

If  the  other  injection  parameters  are  kept  constant  and  u  is  allowed  to  vary,  then 
the  distance  of  closest  approach  will  be  smaller,  the  larger  is  the  wilue  of  u  . 


'  u2 

>/  C)  )2  *  T  -  W  (  Pr 
* 


0  (50) 
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The  trichotomy  in  the  forbidden  zone  types  is  now  shown  to  be  a  reflection 
of  the  trichotomy  in  the  meridional  plane  potential  types.  Three  qualitatively 
different  potentials  exist  corresponding  to  whether  y  is  positive,  negative  or 
zero.  These  cases  will  now  be  discussed,  using  the  paraxial  approximation. 

The  equatorial  plane  traces  of  these  potentials  are  plotted  in  Fig.  1  fora 
representative  value  of  v  .  V  =  0,  in  this  case 

(P,  C)  -  o2/8(1  +  C2)2  (51) 

where  *Q  denotes  the  paraxial  approximation  tow  and  the  superscript  indicates 
the  range  of  values  of  y  .  The  potential  exhibits  a  simple  variation;  it  increases 
as  the  square  of  p  for  constant  C  ,  and  decreases  away  from  the  equatorial  plane 
as  T2  (  O  for  constant  o  .  y  >  0  For  this  case 


+  (0,  5)  =  1/2  [y/p  -  P/2(  1  +  C2) 


(52) 


This  potential  has  a  zero  on  the  line  of  force  given  by 

P  a  (p,  C)  *  y  =  p2/ 2  ( 1  +  r2)  (53) 

Since  this  line  corresponds  to  the  allowed  zone  for  u  =  0,  it  is  seen  that  for 
positive  y  a  line  of  force  is  a  possible  trajectory  in  the  limit  of  infinitely  small 
velocity.  The  centrifugal  potential  contributes  strongly  as  p  -  0,  and  the 
vector  potential  contributes  strongly  as  p  -  ao.  y  <  0,  here 


»Q  (P#  C)  =  1/2 


y/p  +  p/2  ( 1  +  C2 ) 


(54) 


This  potential  form  has  no  zeros,  but  possesses  a  local  minimum  on  the  line  of 
force 

P  o(p,  C)  =  V  -  p2/2(1  +  C2)  (55) 
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because  the  derivative  with  respect  to  p  is  zero  there.  Again  the  contributions 

from  the  centrifugal  and  vector  potentials  show  up  as  in  the  positive  y  case . 

The  difference  between  the  positive  and  negative  potentials  is 
•  + 

*  -  w  -  7/(  1  +C2  )  (56) 

o  o 

_  + 

hence  the  *  surface  is  just  the  w  surface  shifted  vertically  by  an  amount  which 
o  o 

varies  from  y  in  the  equatorial  plane  to  zero  at  r,  -  oo  (see  Fig.  2,  where  this 
difference  has  been  plotted  along  the  line  of  force  given  by  equation  (  55  )  for 
the  same  value  of  V  used  for  the  plots  of  Fig.  1  ).  This  result  is  useful  in  the 
sequel . 

The  forbidden  zone  behavior  is  now  easier  to  understand.  Consider  the 
cases  of  Fig.  1 .  An  allowed  zone  reaches  the  equatorial  plane  in  all  three  cases 
for  values  of  u*  / 2  >  0.0198  =  y ,  while  for  u*  /2  <  0.0198  =  y  only  the 
zero  and  positive  y  zones  penetrate  to  the  equotorial  plane;  this  is  in  accord 
with  equation  (  38  ) . 


D.)  THE  PARAXIAL  TRAJECTORY  EQUATION 


As  a  result  of  the  decoupling  of  the  meridional  and  azimuthal  motion,  it  is 
possible  to  find  directly  a  differential  equation  which  the  meridional  plane  tra¬ 
jectory  satisfies  and  which  includes  only  (  p  ,  £  )  as  variables.  The  Jacobi  form 

o 

of  the  Principle  of  Least  Action  (see  Goldstein  )  here  takes  the  form 


a 

C2 


(h) 


1/2 


1/2 


dC  ■  0 


(57) 
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where  the  following  definition  has  been  introduced 

P  a  u2  -  2w  (  p ,  C  )  (  58  ) 

The  Euler  equotion  for  the  variational  principle  (  57 )  yields  the  following 
differential  equation  for  the  meridional  plane  trajectory: 

p”  =  (  1  +p*  )/2p[0p/dp)  -  P*  (bp  /dp)  ]  (59) 


Where  the  prime  denotes  differentiation  with  respect  to  C  •  Since  this  equation 
ie  too  comp!  icated  to  solve  in  its  full  generality,  the  reason  for  the  use  of  the 
paraxial  approximation  in  studying  the  trajectories  becomes  evident. 

The  paraxial  trajectory  equation  is  obtained  from  equation  (  59)  by 
introducing  the  following  assumptions: 


(I)  Y  =  0 

(II)  C?  <<  1 

(  III  )  pn  «  1 

(IV)  J  >^a2(D,C) 

Thus  the  paraxial  trajectory  equation  for  the  bell-shaped  field  becomes 

p"2  =  p/4u2  (1  +  C2)2 

The  solution  to  this  equation  was  obtained  by  Glaser^®  in  the  form 


=  c/l  +  '2  sin  I  „/l  +  1/4  u2  cot-1 


L-v 


^  +  -f 


(60) 


(61  ) 


where  C  and  q>  are  functions  of  the  injection  parameters.  For  particles  injected 
on  the  axis  with  injection  parameters  (  po  ,  ,  o'  ,  u  )  this  has  been  shown  to 


take  the  form 


11. 


'  J  1  +  r  2/l  +  4u2  'J  1  +  c2  sin  I  \/ 1  +  I, 


P  =  2  >j  p 


Only  this  type  of  trajectory  will  be  considered  here. 


-1 


-1, 


l/4u2  cot  C  -  cot  ^ 
(62) 
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A  few  pertinent  features  of  these  trajectories  will  be  quoted  here  before 
proceeding  to  a  consideration  of  the  injection  problem  ( see  MacLeod ^  ) . 
The  trajectories  are  contained  within  the  envelope  curve 

p,  -  2  up'o  ,yT7^r  *  4u»/ 

The  trajectory  touches  the  axis  N  times  at  the  point 
P  n»2u 


(63) 


m  a  cotivr+4?  +  cot  c 


n  is  given  by 

r  \  ■*. 

INT  \[\  +  l/4u2  )  (  1  -(cot-1  ir  ) 

where  INT  (x)  signifies  the  largest  integer  less  than  x.  The  trajectory  slope 
may  be  shown  to  have  local  maxima  at  the  points  ”  ,  where  the  slope  is 

»■„  -  <-»n+' 


(64) 


(65) 


(66) 


E.)  DEFINITION  OF  INJECTION  PROBLEM 

It  is  now  possible  to  define  the  injection  problem  more  closely.  Since 
for  the  given  type  of  field  in  region  in  and  around  the  equatorial  plane  becomes 
the  high  field  region,  as  the  term  is  understood  here,  the  injection  problem  may 
be  stated  as  follows: 

To  determine  those  values  of  the  injection  parameters  (tQ;  pQ;  90;  rQ;  pQ;  &0;  £0  ) 

which  will  allow  the  particle  to  attain  the  equatorial  plane. 

% 

It  is  clear  that  the  initial  time  may  be  chosen  arbitrarily  ( the  field  is 
magnetostatic  and  no  particle-particle  phase  effects  are  considered  here  ),  thus 
may  be  taken  to  be  zero . 
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Furthermore,  the  field  it  axially  symmetric,  hence  the  value  of  @0 

has  no  effect  on  the  trajectory  behavior  in  the  meridional  plane  and  may 

be  set  equal  to  zero  with  no  loss  of  generality. 

In  addition  to  the  foregoing  simplifications,  it  will  prove  convenient 

to  introduce  the  trio  (  u,  y  /  t )  Instead  of  (  p  ,  0O ,  £  ),  where 

o  o 

♦  =  tan  1  (dp/drc)  a  tan"1  ( P*0  /  )  (67) 

Thus  the  aim  is  to  find  acceptable  values  of  the  set 

(  P-  #  CQ,  u  ,  Y,  t)  (68) 

o  o 
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IV.)  RESULTS 

It  It  convenient  to  dlsttnguieh  two  eatet:  zero  and  non -zero  y  . 
Qualitative  result!  on  the  specification  of  the  injection  parameters  will 
be  given  for  both  cases  through  a  discussion  of  the  relevant  forbidden  zones, 
and  quantitative  results  will  be  given  in  addition  for  the  zero  v  case  by 
utilizing  the  paraxial  trajectories. 

Case  i:  y  *  0 

These  forbidden  zones  depend  only  on  the  parameter  u  .  Several 

of  the  paraxial  zones  are  depicted  in  Fig .  3.  The  allowed  region  for 

a  particle  of  parameter  u  extends  from  the  axis  to  the  appropriate  u 

curve.  The  axis  is  on  allowed  region  for  all  values  of  u  for  particles 

injected  with  p  =  0,  t  =  0,  w.  However,  as  soon  as  becomes 
o 

non-zero  the  value  of must  be  reckoned  with . 

For  small  values,  the  chances  of  a  remotely  injected  particle 
(  C  2  )  penetrating  to  the  equatorial  plane  are  poor  unless  t  is 
quite  small,  although  there  is  a  trajectory  compression  mechanism  at 
work  which  increases  the  chances  over  the  purely  geometrical  values. 

As  u  increases  the  acceptance  angle  (  defined  to  be  the  maximum  value 
of  i  whichallows  equatorial  plane  attainment  )  at  any  given  C 

o 

increases.  It  is  not  possible,  however,  to  increase  u  indefinitely 

and  keep  increasing  the  acceptance  angle.  According  to  the  paraxial  zones 

the  equatorial  plane  zone  radius  increases  indefinitely  as  u  is  increased; 
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this  it  an  error  due  to  the  approximation  used.  Although  the  paraxial  zones 
do  not  show  it,  the  field  source  is  located  at  p  *  1 .0  in  the  equatorial 
plane,  hence  the  equatorial  zone  radius  should  approach  1 .0  asymptotically 
as  u  -*  ao;  only  for  infinite  velocity  can  the  particle  reach  the  field  source. 

Thus  the  error  in  the  zone  boundaries  introduced  by  using  the  paraxial  ap- 
p»0K (motion becomes  more  and  more  serious  as  u  increases.  A  consideration 
of  the  second  term  in  the  vector  potential  expansion,  equation  (  25  ),  shows 
that  this  term  contributes  10%  to  the  vector  potential  at  p  «  0.63  (  in  the 
equatorial  plane  ).  The  toleranco  of  the  error  in  the  specification  of  the 
zone  boundary  position  varies  with  the  application,  but  a  10%  error  will 
be  taken  here  to  be  the  maximum  allowable  value.  In  effect,  this  limits 
the  range  of  u  which  may  be  used  with  reasonable  confidence  becouse  any 
zone  which  has  an  equatorial  radius  o  >  0 .63  has  at  least  a  10%  error  as¬ 
sociated  with  it.  Since  the  equatorial  radius  is  given  by 

p  eq  =  2  u  ( 69 ) 

u  is  limited  by  the  relation 

u  £  P  eq/ 2  es  0.30  (10%  error)  (70) 

It  should  be  noted,  though,  that  this  limitation  applies  strictly  only  to  those 
particles  which  approach  the  zone  boundary  in  or  near  the  equatorial  plane. 

The  choice  of  the  romaining  injection  parameters  (pQ,  r0  /  '))  can  greatly 
change  the  limit  (  70  ) .  The  extreme  example  of  this  is  axial  injection  with 
♦  *0  or  it;  any  value  of  u  is  acceptable . 
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The  ute  of  the  paraxial  trajectories  is  limited  in  like  manner  by  the 
assumptions  (  I  )  -  (  IV )  above .  Assumptions  ( ill )  and  (  IV  )  are  the 
most  restrictive  of  the  group.  According  to  (  IV )  the  trajectories  must 
not  approach  the  forbidden  zones  closely.  At  the  very  least  this  means 
that  the  trajectory  envelope  must  not  intersect  the  zone  boundary.  This 
gives  rise  to  a  condition  on  the  injection  parameters.  Since  the  envelope 
is  similar  to  the  zone  boundary  but  of  smaller  curvature  everywhere  (  compare 
equations  ( 63  )  and  (  35  ) ),  if  the  equatorial  envelope  radius  is  smaller 
than  the  equatorial  zone  radius  the  curves  will  not  cross  anywhere .  Condition 
(  IV )  thus  becomes  in  effect 

Pe  <<  PQ  (  in  the  equatorial  plane )  (71) 

or 

P*  V  1  +  r}  /\  +  4u*  «  1  (72) 

o*  o 

which  is  satisfied  by  small  p^  ,  small  CQ  or  large  u  •  It  is  interesting  that 
this  latter  form  of  condition  (  IV  )  turns  out  to  be  similar  to  condition  (  III  ). 

Since  the  maximum  value  of  pQ  is  at  the  nth  axis  touchpoint,  according  to 
equation  (  66  ),  condition  (  III )  may  be  formulated  as  follows: 

"oV  I  +  ?„/'*  t’n  «  '  (73) 

which  is  satisfied  by  small  pQ  ,  sma  1 1  C  Q  or  la  rge  Cn  *  The  worst  case 
corresponds  to  the  particle  having  one  of  its  axis  touchpoints  near  the  field 
center;  then 

Cn  -  0  (74) 
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and 


P* 

o 


V' 


l+c2  «  1 

O 

hence  condition  (  III  )  is  generally  more  stringent  than  condition  (  IV  ), 
although  for  small  values  of  u  the  difference  becomes  si  ight. 

Although  it  may  seem  that  condition  (  II )  is  the  most  stringent  of  all, 
this  is  not  the  case .  The  purpose  of  imposing  (  It )  was  solely  to  ensure 
rapid  convergence  of  the  vector  potential  series;  the  series  may  converge 
rapidly  without  {  II  ) .  Such  is  the  case  for  certain  regions  of  the  meridi¬ 
onal  plane,  in  fact  for  the  region  within  which  the  contribution  from  the 
second  term  is  less  than  a  specified  small  fraction  of  the  leading  term, 
say  0.10.  The  boundary  of  this  region  may  be  obtained  by  equating 
the  second  term  to  0.10  times  the  leading  term  and  thereby  reduces  to 
the  equation 

0  =  0.632 


H  CJ 


7||  -3 

for  C  >>  1  this  becomes 

o  =  0.365  £ 

a  considerable  weakening  of  condition  (  II  ),  and  the  condition  that  pe 

be  less  that  this  becomes  _ _ ^ 

I 

,2 


0  P 


O  1  +  cy  1  +4uz  ^  0.365 


(75) 


(76) 


(77) 


(78) 
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If  u  is  small,  this  is  easily  satisfied;  if  u  is  large,  (  78  )  takas  the  form 

>0  (  '  +  Co  )  ?  °-730 

which  is  less  restrictive  than  equation  ( 75  ) . 

Thus  the  conditions  (  II )  -  (  IV  )  will  a  I*  be  satisfied  if  equation 

(75  )  is  valid;  this  is  the  basic  restriction,  and  it  is  serious,  though  not 

of  unexpected  character.  For  remote  injection  (  75  )  becomes 

o’  «  l/C 

o  o 

i.e., 

tan  *  ^  \/r 

o 

or  roughly, 


(79  ) 


(80) 

(81) 


*  «  l/CQ  (82) 

The  injection  angle  must  be  considerably  less  than  the  angle  subtended 
by  the  inner  edge  of  the  field  source  (  p  -1.0,  £  =  0  );  this  was 

to  be  expected . 

Thus  the  injection  parameter  ranges  for  axial  injection  are  os  follows 
(  with  the  above  qualifications  ): 

u  <  0.30  (  83  ) 

tan  ♦  «  \/Qq  (84) 

and  r ^  is  arbitrary,  but  its  choice  is  influenced  by  condition  (  84  ) . 
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Cosa  II;  y  f1  0 

The  basic  question  to  be  answered  here  is  whether  there  seems  to  be  any 
advantage  to  non-zero  y  injection  which  would  suggest  that  it  be  used  in 
spite  of  the  less  detailed  knowledge  available  about  trajectory  properties. 

The  answer  to  this  question  seems  to  be  affirmative. 

First  of  all  there  is  a  qualitative  difference  between  positive  and 
negative  v  injection  (  see  Figs.  4  and  5  for  representative  zones  ) .  The 
allowed  zones  for  positive  y  extend  between  the  two  branches  of  the  zone 
boundary  curves  and  include  the  line  of  force.  The  allowed  zones  for 
negative  y  are  to  the  right  of  the  zone  boundary  curvos  and  include  only 
part  of  the  line  of  force  (  no  two-branch  zones  are  plotted  for  negative  y  , 
although  they  exist  when  condition  (  85  )  below  is  fulfilled  ).  No  particle 
of  negative  y  may  reach  the  equatorial  plane  regardless  of  the  other  injection 
parameter  values  unless 

u2  >  2y  (85) 

This  is  emphasized  by  the  plot  of  v  -  along  the  line  of  force  which  is  given  in 
Fig.  2.  Asa  particle  moves  along  the  line  of  force  toward  the  equatorial 
plane  i  -  increases  rapidly  and  eventually  forces  the  particle  to  retreat 
unless  condition  (  85  )  is  satisfied.  Even  when  this  condition  is  satisfied,  part 
of  the  particle  kinetic  energy  is  expended  in  overcoming  this  equatorial 
barrier  potential .  On  the  other  hand,  if  the  other  injection  parameters  have 
appropriate  values,  a  particle  of  positive  y  is  highly  preferable  to  Injection 
with  negative  y  . 
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Secondly,  the  forbidden  zone  is  approached  as  the  azimuthal  velocity 
is  increased  when  the  particle  nears  the  equatorial  plane.  Therefore,  it  would 
seem  advantageous  to  inject  with  zero  angular  velocity  in  order  to  increase  the 
chances  for  the  particle  to  penetrate  to  the  equatorial  plane. 

Thirdly,  injection  along  a  line  of  force  would  seem  to  offer  the  best 
chance  for  equatorial  penetration  because  the  meridional  plane  force  on  the 
particle  is  zero  on  the  line  of  force .  Thus  the  following  condition  must  be 
satisfied  for  line  of  force  injection: 

P*/2(1  +  C02)  (86) 

lastly,  it  is  desirable  to  use  a  relatively  small  value  of  u  to  allow 
the  magnetic  field  to  influence  the  particle  in  the  high  field  region.  A  bal¬ 
ance  must  be  struck  between  too  large  a  value  of  u  (  too  little  field  influence  ) 
and  too  small  a  value  of  >  (  too  much  field  influence  ) . 
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